Određivanje uopštene kompleksne funkcije
assume(0 < R & 0 < C & 0 < L )
zamena = L == C*R^2
zamena = 
vrednosti = [C == 1, R == 1]
vrednosti = 
H(s) = simplify(subs(R/(R + s*L + 1/(s*C) + R), lhs(zamena), rhs(zamena)))
H(s) =

Amplitudska i fazna karakteristika
Hjw(w) = simplify(subs(H(s), s, 1i*w))
Hjw(w) =

Aw(w) = simplify(abs(Hjw(w)))
Aw(w) =

phi(w) = angle(Hjw(w))
phi(w) =

Crtanje amplitudske i fazne karakteristike
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
fplot(w, subs(phi(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
Određivanje propusnog opsega
Aw_num(w) = subs(Aw, lhs(vrednosti), rhs(vrednosti))
Aw_num(w) =

wg = solve(diff(Aw_num,w) == sym(0), w, 'Real', true)
Aref = Aw_num(wg)
Aref =

w3dB = simplify(solve(Aw(w) == Aref/sqrt(2)))
w3dB =

BandPass3dB = w3dB.'
BandPass3dB =

wg1 = subs(w3dB(1), lhs(vrednosti), rhs(vrednosti))
wg1 = 
wg2 = subs(w3dB(2), lhs(vrednosti), rhs(vrednosti))
wg2 = 
Označavanje propusnog opsega u grafiku
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
rectangle('Position',[double(wg1), 0, double(wg2-wg1), double(Aref/sqrt(2))],...
'FaceColor',[1 0.9 0.9],...